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We put forward a strategy to encode a quantum operation into the unmodulated dynamics of a quantum net¬ 
work without the need of external control pulses, measurements or active feedback. Our optimization scheme, 
inspired by supervised machine learning, consists in engineering the pairwise couplings between the network 
qubits so that the target quantum operation is encoded in the natural reduced dynamics of a network section. 

The efficacy of the proposed scheme is demonstrated by the finding of uncontrolled four-qubit networks that 
implement either the Toffoli gate, the Fredkin gate, or remote logic operations. The proposed Toffoli gate is 
stable against imperfections, has a high-fidelity for fault tolerant quantum computation, and is fast, being based 
on the non-equilibrium dynamics. 


INTRODUCTION 

Computational devices based on the laws of quantum me¬ 
chanics hold promise to speed up many algorithms known to 
be hard for classical computers m. The implementation of 
a full scale computation with existing technology requires an 
outstanding ability to maintain quantum coherence (i.e. iso¬ 
lation from the environment) without compromising the abil¬ 
ity to control the interactions among the qubits in a scalable 
way. Among the most successful paradigms of quantum com¬ 
putation, there is the “circuit model”, where the algorithm 
is decomposed into an universal set of single- and two-qubit 
gates m, and, to some extent, the so-called adiabatic quan¬ 
tum computation (AQC) a where the output of the algorithm 
is encoded in the ground state of an interacting many-qubit 
Hamiltonian. A different approach ||4] is based on the use of 
always-on interactions, naturally occurring between physical 
qubits, to accomplish the computation. Compared to the cir¬ 
cuit model, this scheme has the advantage of requiring min¬ 
imal external control and avoiding the continuous switch off 
and on of the interactions between all but two qubits; while 
compared to AQC it has the advantage of being faster, being 
based on the non-equilibrium evolution of the system. Quan¬ 
tum computation with always-on interactions is accomplished 
by combining the natural couplings with a moderate exter¬ 
nal control, e.g. with a smooth shifting of Zeeman energies 
0, via feedforward techniques using measurement based 
computation ||2l or quantum control 0|9l. Most of these ap¬ 
proaches are based on the assumption that the natural cou¬ 
plings are fixed by nature and not tunable, while local inter¬ 
actions can be modulated with external fields. However, the 
amount of external control required can be minimized if the 
couplings between the qubits can be statically tuned ifTOl . e.g. 
during the creation of the quantum device. 

The recent advances in the fabrication of superconducting 
quantum devices has opened up to the realization of inter¬ 
acting quantum networks. In a superconducting device, the 
qubits are built with a Josephson tunnel element, an induc¬ 
tance and a capacitor CD, while local operations and mea¬ 
surements are performed by coupling the qubit to a resonator 
CD- The interactions can be designed using lithographic 


techniques by jointly coupling two qubits via a capacitor d 
or an inductance oa, and can be modeled via an effective 
two-body Hamiltonian 'Y^^JaO'a®cra ifTSlfTbl . where ctq are 
the Pauli matrices. Because of the flexibility in wiring the 
pairwise interactions among the qubits, it is possible to ar¬ 
range them in a planar graph structure, namely a collection 
of vertices and links, where the vertices correspond to the 
qubits and the links correspond to the 2-body interactions be¬ 
tween them. Moreover, thanks to the development of three- 
dimensional superconducting circuits ini, it may be possi¬ 
ble in the near future to wire also non-planar configurations, 
namely a general qubit network. 

Motivated by the above arguments, we consider the ques¬ 
tion whether it is possible to encode a quantum algorithm into 
the unmodulated dynamics of a suitably large quantum net¬ 
work of pairwise interacting qubits. This would be extremely 
interesting, as it would enable quantum computation by sim¬ 
ply “waiting”, without the need of continuously applying ex¬ 
ternal control pulses or measurements. Even when sequen¬ 
tial operations cannot be avoided, our scheme can enable the 
in-hardware implementation of recurring multi-qubit opera¬ 
tions of a quantum algorithm (see e.g. Fig{T]), such as quan¬ 
tum arithmetic operations QSl, and possibly also the quan¬ 
tum Fourier transform or error correcting codes HI. We fo¬ 
cus on two-body interactions, since they are the most com¬ 
mon in physical setups, and we consider an enlarged net¬ 
work where auxiliary qubits enrich the quantum dynamics. 
The important question analyzed in this paper is the follow¬ 
ing: given a target unitary operation Uq on a given set of 
qubits Q, we consider an extended network QuA where A 
is a set of auxiliary qubits {ancillae), and we ask whether it is 
possible to engineer the pairwise interactions in QUA, mod¬ 
eled by the time-independent Hamiltonian Hqa, such that 
e^^^Q^=UQ(S)VA after some time t {Va may be an extra uni¬ 
tary operation on the auxiliary space). More generally the tar¬ 
get operation can depend also on the ancillae initial state: if 
gjir/QA— C/q”^( 8)|A„)(A„|, where {|A„)} form a basis of 

the ancillae Hilbert space and, e.g., Uq=Uq \ then the target 
operation is implemented when A is initialized in |Ai). Our 
method is particularly useful for implementing quantum gates 
which requires /c-body interactions (fc>2), such as the Toffoli 
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FIG. 1. Schematic transposition of a quantum circuit to a trained quantum network, (a) An example 7-qubit circuit where the gates Gi and 
G 2 are sequentially applied. Many important circuits can be cast into the scheme (a), such as those for quantum arithmetics (H. (b) Quantum 
network implementation of (a): each qubit network in the green boxes implements either Gi or G 2 on the input/output qubits (the three qubits 
in the bottom row). The quantum bus sequentially transfers the state of three qubits from the register, upon which the gate Gi or G 2 has to 
be applied, to the input/output qubits of the gate network. After the gate network has implemented its transformations, the state is transferred 
back to the original three-qubits of the register. 


or Fredkin gates |[Tl[T9l|20l where for any 2-local 

Hq, and for remote logic, namely for applying a gate to qubits 
which are not directly connected but are rather interacting via 
intermediate systems. Our approach is completely different 
from the simulation of fc-local Hamiltonians with pairwise 
interactions discussed in the AQC literature ll^ 1^ . being 
based on the unmodulated dynamics. Moreover, being based 
on unmodulated (time-independent) interactions and ancillary 
qubits, it is significantly different from quantum optimal con¬ 
trol na. 


Our quantum network design procedure is inspired by su¬ 
pervised learning in feedforward networks EH, where the 
training procedure involves the optimization of the network 
couplings (i.e. the weights between different nodes) such 
that the output corresponding to some input data has a de¬ 
sired functional form (e.g. for data classification). Although 
there are many recent developments about using a quantum 
device to speed-up machine learning algorithms ll25H29l or 
storing data ll^ . our optimization procedure is entirely clas¬ 
sical, but specifically developed for quantum hardware design. 
Our scheme is completely different from other recent propos¬ 
als EiHia because it avoids measurements or active feed¬ 
backs and requires minimal external control. 


RESULTS 


Supervised quantum network design 


Supervised learning is all about function approximation: 
given a training set {(/i, Oi), (/ 2 , O 2 ),. •. }, namely a col¬ 
lection of inputs Jfc and the corresponding known outputs Ok, 
the goal is to find a function / with two desired properties: i) 
Ok—f{Ik) for training pair; ii) / should be able to infer 
the unknown output of an input not contained in the training 
set. In classical feedforward networks, the function / is ap¬ 
proximated with a directed graph organized in layers, where 
the first layer is the input register and the last one encodes the 
output. The value s). of the fc-th node in layer £ is updated 
via the equation ’ where is 

an appropriate (typically non-linear) activation function and 


A 


(e-i) 

kj 


is the weight between node k in layer £ and node j 


in £—1. The training procedure consists in finding the opti¬ 
mal weights A by minimizing a suitable cost function such as 

c^EklOk-fih)?- 

A quantum network consists on the other hand of an undi¬ 
rected graph (y, E) of vertices V and links E described by a 
2-local Hamiltonian 


{n,m)£E 






nGV a 


where cr“, a=x,y,z, are the Pauli matrices acting on qubit 
n and, to simplify the notation, we call 
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the set of parameters. The vertices are composed of 
two disjoints sets V=QUA where Q consists of regis¬ 
ter qubits and A of auxiliary qubits. Given a sep¬ 
arable initial state |^q) 0 the time evolution ac¬ 

cording to Hamiltonian Q generates a quantum channel 

Cl ^a[|'0q) (V'Q|]=TrA[e"*’^* \tpQ) (V'qI® \iPa) 

on subsystem Q - since we are interested in a fixed opera¬ 
tional time t for simplicity we set t=l, reabsorbing t into the 
definition of the definition of H. Depending on the flexibility 
of the experimental apparatus in reliably initializing the aux¬ 
iliary qubits, one can add {iPa) to the set A. Network design 
consists in the following procedure: given a target unitary op¬ 
eration Uq that we want to implement, the goal is to find the 
parameters A, if they exist, such that £x[pq]=UqPqUq for 
any pQ. To simplify the notation we assume that the gate out¬ 
put is encoded in Q but it is straightforward to generalize the 
formalism when the output sites differ from the input ones. 

Motivated by the similarity with classical supervised learn¬ 
ing, where the weights A are tuned to maximize the ability 
of the network to reproduce a known output given the corre¬ 
sponding input, we create a training set T with a random set of 
initial input states. For each input Ipj) sT the expected known 
output is Uq \tp), while the output of the network evolution is 
£’x[|'*A) i'fpW- The “learning” procedure involves the minimiza¬ 
tion of the difference between the output of the network and 
the expected output, and corresponds to the maximization of 
the fidelity 


■^= E j^p = {mh£xmmuQ\p^}. (2) 

Wer'' ' 

If the average is performed over all possible states, then Eq. © 
can be substituted by the average gate fidelity J- = J Tjjjdip 
where the formal integration can be explicitly evaluated cni 
|3]|35l yielding 


T = 


1 


D + 1 






(3) 


where £"^'’"^ = {qi\£x[\qk){qi\]\qj), Uij={qi\UQ\qj) and 

form the computational basis of the £)-dimensional 
Hilbert space of qubits Q. The typical value of the fidelity 
for a random non-optimal evolution of the qubit network is 
F=D~^, obtained using Haar integration techniques 1^ . 
This value is independent on the details of the ancillae, since 
it depends only on the dimension of the target Hilbert space, 
and provides an estimate for the initial fidelity of an untrained 
network. 

The gate learning procedure corresponds to a global max¬ 
imization of the fidelity ([^. However, because of the many 
parameters in the Hamiltonian Q, IF can have many local 
maxima making the global optimization extremely compli¬ 
cated. Since most global optimization algorithms introduce 
stochastic strategies, rather than introducing unphysical ran¬ 
dom jumps, we take advantage of the explicit stochastic nature 
of the problem (IF is a uniform average over random states) 
and we propose the following learning algorithm to design the 
interactions of the quantum network. 



FIG. 2. Average fidelity F and fidelity for some random states tp 
for implementing the Toffoli gate in a 4 qubit network (see discussion 
in the text). All the parameters are set to the optimal ones, except J-j§ 
in the abscissa. The region around the only global peak is filled in 
green. 


1: Choose an initial parameter set A (e.g. at random), and 
choose an initial learning rate e; 

2: repeat 

3: generate a random \ip)', 

4: update £ times the coupling strengths as 

\^X + eVx{mh£xmmUQ\i^); (4) 

5: decrease e (see Methods); 

6: until convergence (or maximum number of operations). 

Specifically, we combine the above algorithm with the max¬ 
imization of the average fidelity (see below) and we observe 
a drastic speedup of the optimization process. The param¬ 
eter £ tunes the number of deterministic steps in the learn¬ 
ing procedure, and can be set to the minimum value £=1, so 
that after each interaction the state is changed, or to a higher 
value. In our simulations we use £=1, for simplicity. Our 
algorithm is an application of the stochastic gradient descent 
(SGD) method ll37l to the maximization of the function (|^. 
In classical feedforward networks, SGD is the de facto stan¬ 
dard algorithm for network training ll24llJ7l and is specifically 
used for large training datasets, when the evaluation of the cost 
function and its gradient are computationally intensive. On 
the other hand, the average in Eq.(|^ can be evaluated explic¬ 
itly over a uniform distribution of an infinite number of initial 
states, giving Eq. Although is easier than F to com¬ 
pute, the major advantage of SGD for quantum network de¬ 
sign comes from its ability to escape local maxima. The cru¬ 
cial observation to show the latter point is that the statistical 
variance over random states NwiF=F’^—F'^ vanishes when 
F=\ (see e.g. EiEa) - indeed, intuitively, since both F and 
F^ are bounded in [0, 1], F can achieve its maximum only if 
£"^=1 for all the states, apart from a set of measure zero. On 
the other hand, if 0<£<1, then Var£>0 and the fluctuations 
can be so high that a local maximum of F may not correspond 
to a maximum of F^ for some state ip. This is indeed shown 
in Eig|^ with a real example for the implementation of the 
Toffoli gate (see the application section below). In Eigj^the 
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average fidelity T has three local maxima at (fc = 1,2,3) 
and a single global maximum at A®^', namely the optimal pa¬ 
rameters, while the fidelities for different random states tjj 
have a more complicated behavior. In view of the argument 
discussed above, all the state fidelities have a global max¬ 
imum at A®^ while, remarkably, at least one fidelity has 
no local maximum at a 5°‘^' . Our stochastic learning algorithm 
uses a gradient descent technique for locally maximising the 
function T^{\). Therefore, if we are around the slopes of a 
local maximum of (say from the previous exam¬ 

ple) and the state lip) is randomly changed to \(p), that local 
maximum may disappear from allowing the algorithm 

to escape from this non-optimal region when the parameters 
are updated via Eq. 0- On the other hand, when the algo¬ 
rithm is probing the neighborhood of a true optimal point for 
which ,F(A)=1, (e.g. A®*' in the previous example), then the 
maximum of does not disappear when the state ip is 

changed, allowing the “climbing” procedure to continue. 

The above stochastic algorithm may be combined with a 
deterministic maximization of Eq.Q. In our simulations we 
use stochastic learning for the initial global span of the param¬ 
eter manifold and, if it reaches a suitably high fidelity (e.g. 
,F>95%), then it is reasonable to suppose that the algorithm 
has found a global maximum. Starting from this point we per¬ 
form a local maximization of Eq. 0 and, if is reached, 
the learning has been successful. Otherwise we repeat the pro¬ 
cedure. 

It is worth emphasizing that given a target gate U, it is 
an open question to understand a priori whether a solution 
may exist for a graph with a certain set of interactions (e.g. 
Heisemberg, Ising, etc.). Unlike in quantum control, where 
given a time dependent Hamiltonian 'H(t) = T-Lq + P{t)V one 
can check in advance whether U = T[exp(—z 'H{t)dt)\ 
for some control profile f3{t): such profile can exist only if 
U is contained in the group associated to the algebra gener¬ 
ated by the repeated commutators of Hq and V. Although no 
complete algebraic characterization is known for our case (see 
however the Methods for a necessary condition) and we have 
to study each problem numerically, in the next sections we 
find some structures which enable the implementation of im¬ 
portant quantum gates. All numerical simulations have been 
obtained in a laptop computer using QuTiP ll38l . 


Application: Toffoli gate 

The Toffoli gate is a key component for many important 
quantum algorithms, notably the Shor algorithm lf3^ , quan¬ 
tum error correction ||20l, fault-tolerant computation iQi, 
quantum arithmetic operations m and, together with the 
Hadamard gate, is universal for quantum computation ED. 
Experimental implementations of this gate has been obtained 
with trapped ions ll42l . superconducting circuits im in, 
or photonic architectures IBl. Toffoli gate is a controlled- 
controlled-not (CCNOT) operation acting on three qubits. It 
can be implemented in a circuit using five two-qubit gates fD, 
or can be obtained in coupled systems via quantum control 



FIG. 3. Network implementing the Toffoli gate. The gate acts on the 
three external qubits (the top ones being the control qubits, and the 
bottom one being the target), and has an additional auxiliary qubit in 
the center. 


techniques EIlISl. Efficient schemes require higher dimen¬ 
sional system (i.e. qudits) EH. On the other hand, the di¬ 
rect implementation using natural interactions is complicated, 
since the Hamiltonian "Hccnot corresponding to the gate, i.e. 
CCNOT=e*^°'="°'^, has three-body interactions which un¬ 
likely appear in nature. 

By applying our quantum hardware design procedure, we 
show that the Toffoli gate can be implemented in a four qubit 
network using only pairwise interactions and constant con¬ 
trol fields. Our findings enable the construction of a device 
which implements the Toffoli gate with a fidelity ,F=99.98% 
by simply “waiting” for the natural dynamics to occur, with¬ 
out the need of external control pulses. We consider a four 
qubit network as displayed in Eig. where the control qubits 
are labeled by the indices 1,2, the target is qubit 3 and the an- 
cilla is qubit 4. We start our analysis by considering a fully- 
connected graph where each qubit interacts with the others 
using XX- and ZZ-type pairwise interactions, as this kind of 
interaction can be obtained in superconducting circuits M- 
Because of the symmetries of the Toffoli gate (see Methods), 
we consider the two control qubits to be equally coupled to 
the target and the ancilla: w=3,4 and simi¬ 
larly we set ■ Moreover, since the Toffoli gate is real, 

we only consider local fields in the X and Z directions and set 
\ip a) = cos Tj It) sinp |t). By combining SGD with the 

maximization of Eq.([^ we find the following optimal param¬ 
eters. 


= -8.940, 
hi = -2.428, 
hi = -19.08, 

p = 0.8182, 


Jll = -4.957, 
hi = Jit, 
hi = -4.267, 
C = 0.0587, 


= -5.657, 
hi = -0.165, 
Jgt = 15.06, 


where the other XX- and ZZ-type interactions not displayed 
in 0 are found to be zero by the learning algorithm, so the 
optimal configuration is the one summarized in Eig|^ where 
the XX coupling is only between qubits 3 and 4. In more 
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physical terms, if the maximal allowed coupling is fixed to 
J/27r«40MHz, then we find a gate time of 60ns and 


= -149.2 MHz, 
= -94.39 MHz, 
= -318.4 MHz, 
hi = -40.52 MHz, 
K = J!z- 


= -82.71MHz, 
73 ^ 4 "^ = 251.3 MHz, 
hi = -71.2 MHz, 
hi = -2.751MHz, 


With the optimal parameters of Eqs. 0,0 we obtain an 
average gate fidelity of 99.98%, above the threshold for 
topological fault tolerance for single- and two-qubit gates, 
while by avoiding the extra phase fixing ^=0 we still obtain 
.F=99.92%. Moreover, our gate fidelity is above the Toffoli 
gate accuracy threshold (755/756~99.87%) for fault-tolerant 
computation in the limit in which Clifford gate errors are neg¬ 
ligible BtI . 

The optimal parameters 0,0 are stable against an imper¬ 
fect tuning of the interactions. Indeed, we considered a pertur¬ 
bation Afe—^Afe-l-erfc, rfcS[ 0 , 1 ] being a random number and e 
being the strength of the static perturbation, and found that 
^■>99.9% if e<0.04 (e<0.7MHz) and .F>99% if e<0.18 
(e<3MHz). 


Application: Fredkin gate 

Fredkin gate is a controlled-swap (CSWAP) opera¬ 
tion acting on three qubits which is universal for re¬ 
versible computation m. We found that this gate can 
be obtained with perfect fidelity (up to the numerical 
precision) in a four qubit network with Hamiltonian Q 
where Jf^=Jf^=13.60 (227.0 MHz), J^3“=-4.712 

(-78.62 MHz), =8.400 (140.2 MHz), 

Jf|=Jf|=11.15 (186.1MHz), /if=1.025 (17.11 MHz), 
hf=TT (54.42 MHz). The values in MHz correspond to a gate 
time of 60ns. Moreover, e“*^=CSWAPi 23 ( 8 )H 4 so the gate 
is independent on the initial state of the ancilla. As for the 
Toffoli gate, this optimal configuration has been obtained by 
starting the training procedure with a fully connected graph 
with all the interactions, so the fact that some interactions are 
zero is a result of the optimization process. 


Application: remote logic 

We study a qubit network which implements a maximally 
entangling gate between two sites which are not directly cou¬ 
pled. Remote logic has been studied extensively in spin chains 
for achieving entangling operations between the boundary 
sites naiioiiiHiiisi, and it is a building block for a proposed 
architecture for solid-state quantum computation at room tem¬ 
perature m- For simplicity we consider a SU(2) invariant 
four qubit network, interacting with a Heisenberg Hamilto¬ 
nian iT= 'Ea=x,y,z where there is no di¬ 

rect coupling between qubits 1 and 4 (Ji 4 = 0 ). Applying our 
learning algorithm we found that the \/SWAP gate, which is 


universal for quantum computation when paired with single 
qubit operations IT], can be achieved between qubits 1 and 4 
with unit fidelity with different choices of Ji 2 =J 24 , 
and J 23 when the initial state of ancillae is (|tT) ~ ITt))/'\/2- 
Given this simplification one can then find a solution analyti¬ 
cally: Ji 2 =a+'nyj ( 2 n)^ — 1 /-\/ 8 , Ji 3 =Q!—tta/( 2 n )2 — 1 /-\/ 8 , 
^ 23 = 0 ! -f (—l)" 7 r, where n is an integer. We find analyti¬ 
cally that irrespective of a the above choice gives perfect fi¬ 
delity. Our strategy has not found any three-qubit configura¬ 
tions which implement a remote \/SWAP gate, so the four 
qubit network is the minimal non-trivial example. Remark¬ 
ably, some of our 4-qubit configuration are more stable to 
noise than the direct implementation of the gate in a two-qubit 
system (namely when Ji 4 = 7 r /2 and the other couplings are 
zero). For instance, if j^optimai 

optimal value for implementing the gate, we found that when 
e is randomly distributed in [0,1/2], then the 4-qubit system 
with n = 1 still has, on average, F ~ 99.1%, while the direct 
two-qubit case has F ~ 98.8%. 


Towards a scalable architecture for quantum computation 

Current architectures for quantum computation, e.g. with 
superconducting qubits m or ion traps ll52l . are based on an 
arrays of interacting qubits which are continuously controlled 
via external pulses to implement the desired operation. This 
approach may suffer from scalability issues because, even as¬ 
suming the ability to maintain quantum coherence for a long 
time, extremely large (classical) control units will be neces¬ 
sary to generate the sophisticated pulse sequences required 
to implement a full-scale quantum algorithm. On the other 
hand, the approach that we have in mind shares more similar¬ 
ities with integrated circuits in nowadays electronics, where 
a set of special-purpose logic units (modules) are wired to¬ 
gether to achieve computation (or other tasks). In our vision, 
different modules can be fabricated with qubit networks de¬ 
signed to produce a specific logic task, namely a quantum 
gate, automatically without the need of external control. As in 
Figj^the different logic and memory units can be reciprocally 
connected using a quantum bus, whose purpose is to transfer 
the qubit states between the quantum registers/memory and 
the input/output qubits of the modules. In Figj^for simplic¬ 
ity the input/output ports of the modules are designed in the 
same physical qubits, although this can be easily extended to 
more general cases. The quantum bus can be realized with dif¬ 
ferent technologies, e.g. with microwave resonators ll53]l . or 
can also be implemented via quantum state transfer in a qubit 
network m. The modules shown in Fig[T] can be designed 
to produce either simple basic operations, like the CNOT or 
the Toffoli gate, or, in principle, they can directly implement 
larger components of a quantum algorithm like the Quantum 
Fourier Transform or error correcting codes IT]. In this re¬ 
spect, to treat systems with many parameters one can easily 
combine our optimization strategy, based on fidelity statistics, 
with metaheuristic strategies 1551 which simultaneously deals 
with many candidate solutions and are known to be fast in 
global optimization with high-dimensional parameter spaces. 
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Moreover, highly optimized deep learning algorithms are al¬ 
ready used to train neural networks with 60 millions of pa¬ 
rameters ||5^ . However, given the difficulty in numerically 
simulating large quantum systems, this approach may be rea¬ 
sonable for networks up to, say, 20-30 qubits. 


DISCUSSION 

Inspired by classical supervised learning, we have proposed 
an optimization scheme to encode a quantum operation into 
the unmodulated dynamics of a qubit register, which is part of 
a bigger network of pairwise interacting qubits. Our strategy 
is based on the static engineering of the pairwise couplings, 
and enables the creation of a quantum device which imple¬ 
ments the desired operation by simply waiting for the natural 
dynamics to occur, without the need of external control pulses. 
Our findings show that machine learning inspired techniques 
can be combined with quantum mechanics not only for data 
classification speed-up ll25l [26l or quantum black-box certifi¬ 
cation El ED, but also for quantum hardware design. 

This paper opens up the topic of encoding quantum gates 
and operations into the unmodulated dynamics of qubit net¬ 
works. Although we have focused on small systems, larger 
networks can be considered using more efficient training 
schemes. These would enable the simulation of larger compo¬ 
nents of a quantum algorithm, since different multi-qubit gates 
can be combined into a unique quantum operation which can 
be simulated in a large quantum network. Moreover, when 
combined with a quantum bus as in Fig. our strategy can 
provide an alternative approach to universal quantum compu¬ 
tation which avoids the decomposition of the algorithm into 
one- and two-qubit gates. Note that most quantum algorithms 
take classical inputs so the extra control required for initializa¬ 
tion demands the further ability to fully polarise globally the 
spins. The latter step is however typically much easier than the 
implementation of entangling gates, which has been consid¬ 
ered in this paper. Moreover, in view of the recent experimen¬ 
tal measures of the average gate fidelity ll59l . it is tempting 
to predict an all-quantum version of our learning procedure 
where IF is not classically simulated, but rather directly mea¬ 
sured. This would require a further highly controlled system 
to infer the optimal parameters of an uncontrolled quantum 
network, which can be used to industrialize the production of 
unmodulated quantum devices implementing the desired al¬ 
gorithm. 

Our results demonstrate the efficacy of the proposed 
scheme in designing four qubit networks which implement the 
Toffoli and Fredkin gates or remote logic operations. The pro¬ 
posed Toffoli gate is fast, has high-fidelity for fault-tolerant 
computation, and only uses static XX- and ZZ-type interac¬ 
tions which can be achieved in superconducting systems M- 
The key advantage of our method is in exploiting all the per¬ 
manent interactions in the qubit network without trying to 
suppress some of them sequentially to implement pairwise 
gates. Moreover, being based on non-equilibrium dynamics, 
our gate is fast; if J/27rR:!40MHz then the total operation time 
is around 60ns which matches the current gate times for single 


and two qubit operations El- 


METHODS 
Learning rate 

The choice of the learning rate e is crucial. If the initial 
learning rate is too small, it might not escape from the differ¬ 
ent “local maximum” points, while if it is too large it will con¬ 
tinue to randomly jump without even seeing the local maxima. 
To maximize the speed and precision of SDG the learning rate 
e has to decrease as a function of the steps, a common choice 
being cocmT^!'^ where m is the step counter El- However, 
when the gradient in Eq.Q cannot be performed analytically, 
one can use more sophisticated techniques where both 
the learning rate and the finite difference approximation of the 
gradient change as a function of m. 


Symmetries 

In the design of the quantum network and its cou¬ 
plings the number of parameters can be drastically re¬ 
duced if the target unitary operation Uq has some 
symmetries, namely if there exists some unitary ma¬ 
trix S such that [[/q,S'q]= 0. This condition requires 
the quantum channel 8x[p]=Tt :to sat¬ 
isfy £a(5'qpS'^)=S'q£’a(p)S'q for each state p, e.g. 
[Tf, 5 'q(8)11a]= 0. Conversely, if the interaction type is fixed 
by nature (for instance, only Ising or Heisenberg interactions 
are allowed), then one has to check whether the Lie algebra 
spanned by the operators in T-L contains the generators of Uq. 


Bottom-up construction: Lie algebraic characterization 

All the numerical results presented in the main text are ob¬ 
tained using a top-down approach; after selecting the interac¬ 
tion types (e.g. XX, ZZ, Heisenberg etc.), the algorithm starts 
with a zero-bias fully connected configuration where all the 
qubit pairs of the network interacts with all possible interac¬ 
tions, each weighted with a different parameter, and different 
local fields. As a result of the training procedure, we found 
numerically that most of these parameters are indeed zero. 
However, for larger networks it is better to use a bottom-up 
approach where one starts with a minimal set of parameters, 
and then adds other parameters until a solution is found. 

To construct a minimal set of parameters one can use a Lie 
algebraic characterization inspired by quantum control. We 
write the Hamiltonian as FL— where Xj are the in¬ 

dependent parameters and Oj the operators. If the parame¬ 
ters are time dependent, then there exist suitable pulses Xj (f) 
such that the dynamics implements the target gate G only if 
log(G) is contained in the algebra generated by the repeated 
commutators [Oj, [Ok, ••■]]■ Since our scheme is based on 
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the particular choice where Xj (t) is constant, the above char¬ 
acterization still provides a necessary condition. As an ex¬ 
ample, we consider the Toffoli gate and the solution Eq. 0 
where C 7 i=CTfcrf, C72=o'f(T|-|-tT|(Tf-|-2cr|, 03=0-^cr|-|-crferf, 
C>4=crfcr4, C>5=(jf-|-cr|, C>7 =ct|, Os=( 7 ^, C79=crf. It is sim¬ 
ple to check that log G (up to an irrelevant constant factor) is 
contained in the algebra generated by the operators Oj, while 
this is not the case if the operator Og is removed from the 
Hamiltonian. Therefore, no solution is possible if A8=0. 

Inspired by the above example the bottom-up approach con¬ 
sists in the following steps; i) based on the symmetries of the 
target gate and on the physically allowed interactions one de¬ 
fines an initial set of operators; ii) other operators are added 


to the set until the dynamical algebra contains log(G'); iii) one 
starts the numerical parameter training to check for conver¬ 
gence (different runs may be required). Until the solution is 
found one then either adds new operators, or change the pre¬ 
vious ones. 
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